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Abstract 

In this paper, we study diagonalizable hyperbolic systems in one space dimension. Based on a 
new gradient entropy estimate, we prove the global existence of a continuous solution, for large 
and nondecreasing initial data. Moreover, we show in particular cases some uniqueness results. 
We also remark that these results cover the case of systems which are hyperbolic but not strictly 
hyperbolic. Physically, this kind of diagonalizable hyperbolic systems appears naturally in the 
modelling of the dynamics of dislocation densities. 
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1 Introduction and main result 
1.1 Setting of the problem 

In this paper we are interested in continuous solutions to hyperbolic systems in dimension 
one. Our work will focus on solution u(t,x) = (u*(t, x))j=i,...,m, where M is an integer, 
of hyperbolic systems which are diagonal, i.e. 

dtu' + a\u)d^u' = on (0, T) x R and for z = 1, ..., M, (P) 
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with the initial data: 

u'{0,x) =ui{x), xeR, for i = 1,...,M. (ID) 

For real numbers a* < P\ let us consider the box 

U = Uf^,[a\n (1-1) 

We consider a given function a = (a*)j=i^...^M : U —>■ M*^, which satisfies the following 
regularity assumption: 



(^1) 



' the function a G C°°(f/), 

there exists Mq > such that for i = 1, ...,M, 
\a'{u)\<Mo for all u E U, 



there exists Mi > such that for z = 1, ..., M, 
laUv) — aUu)\ < Mi\v — u\ for all v,ueU. 



We assume, for all u G M , that the matrix 

(aj(«))i,j=i,...,M, where a[ 
is non-negative in the positive cone, namely 






{H2) 



for all u E U, we have 

^ ^i^ja'j{u) > for every ^ = (^i, ..., 6/) ^ [0, +00^ 

i,j=l,...,M 



M 



In (IIDI) . each component Uq of the initial data Uq = (mq, ■ ■ ■ , Uq^) is assumed satisfy the 
following property: 



ul G L" 



(i/3) < Mq is nondecreasing. 



d^ui G L log L 



for i = 1, 



■,M, 



where LlogL(]R) is the following Zygmund space: 

LlogL(M) = \f EL^iR) such that /" |/| In (1 + |/|) < +00 
This space is equipped by the following norm: 



iiogL(M) = inf |a > : y^ln (l + \S\<i^^ 
This norm is due to Luxemburg (see Adams [H (13), Page 234]). 

Our purpose is to show the existence of a continuous solution, such that u\t, ■) satisfies 
(if 3) for all time. 

1.2 Main result 

It is well-known that for the classical Burgers equation, the solution stays continuous 
when the initial data is Lipschitz-continuous and non-decreasing. We want somehow to 
generalize this result to the case of diagonal hyperbolic systems. 

Theorem 1.1 (Global existence of a nondecreasing solution) 

Assume {HI), {H2) and {H3). Then, for all T > 0, we have: 

i) Existence of a weak solution: 

There exists a function u solution of (f^)-( f/.DI] (in the distributional sense), where 



ue [L°^((0,T) xM)]A^n[C([0,T);LlogL(M))]^ and d^u e [L°°((0,T); LlogL(M))] 



M 



such that for a.e t G [0,T) the function u{t,-) is nondecreasing in x and satisfies the 
following L°° estimate: 

\\u\t, OIIl-cr) < ||mJ,||l-(r), /or z = 1, . . . , M, (1.2) 

and the gradient entropy estimate: 



f Yl f{d.u\t,x))dx+ /V E 



a[j{u)dxu\s,x)dxu\s,x) dx ds < Ci, 

(1.3) 



where 



f(^)- / a;ln(x) + i if x>l/e, . 

^^^'~ \ if 0<a;<l/e, ^^■^> 

and Ci(T, M, Ml, ||-uo||[L°o(R)]Af , ||f^xMo||[LiogL(R)]*0- 

a) Continuity of the solution: 

The solution u constructed in (i) belongs to C([0,T) x M) and there exists a modulus of 
continuity uj[5, h), such that for all (t, x) G (0, T) x M and all 6,h > 0, we have: 

\u(t + 6,x + h) - u(t,x)\ <C2Uj(5,h) with Lu(6,h) = ——. tt + m 7T- (1-5) 

Hs + 1) Hh + 1) 

where C2{T,Mi,Mo, ||mo||[l°°(r)]m, ||<9xMo||[LiogL{R)]M)- 



Remark 1.2 

Here, we can easily extend the solution u of ^)- l^IH) . given by Theorem \l.l\. on the time 
interval [0, +00). 

Our method is based on the following simple remark: if the initial data satisfies (-^3) 
then the solution satisfies (-^3) for all t. What seems new is the gradient entropy in- 
equality. The prove of Theorem 11.11 is rather standard. First we regularize the initial 
data and the system with the addition of a viscosity term, then we show that this regu- 
larized system admits a classical solution for short time. We prove the bounds (11. 2p and 
the fundamental gradient entropy inequality (11.31) which allow to get a solution for all 
time. Finally, these a priori estimates ensure enough compactness to pass to the limit 
when the regularization varnishes and to get the existence of a solution. 



Remark 1.3 

To guarantee the L log L bound on the gradient of the solutions. We assumed in (H2) a 
sign on the left hand side of gradient entropy inequality ^L3^. 



In the case of 2 x 2 strictly hyperbolic systems, which corresponds in ([P]) to the case 
oi a}{u^,u^) < a^{u^,u^). Lax [30] proved the existence of smooth solution of (|P11-( IId11 . 
This result was also proven by Serre [Ml Vol II] in the case of M x M rich hyperbolic 
systems (see also Subsection 11.41 for more related references). Their result is limited 
to the case of strictly hyperbolic systems, here in Theorem 11.11 we treated the case of 
systems which are hyperbolic but not strictly hyperbolic. See the following Remark for 
a quite detailed example. 



Remark 1.4 (Crossing eigenvalues) 

Condition lll.S\) on the eigenvalues is required in our framework (Theorem \l.l\) . Here is 
a simple example of a 2 x 2 hyperbolic but not strictly hyperbolic system. We consider 
solution u = (u^,M^) of 



on(0,r)xM. (1.6) 



dtu^ + cos{u^)dxU^ = 0, 

dfu'^ + u^sin{u'^)dxU^ = 0, 
Assume: 

i) u^{—oo) = 0, u^{+oo) = 1 and d^u^ > 0, 

a) n^(-oo) = -|, M^(+oo) = I and d^u^ > 0. 

Here the eigenvalues \i{u^,u^) = cos{u^) and \2{u^,u'^) = u^sin{u^) cross each other 
at the initial time (and indeed for any time). Nevertheless for a'^{u^,v?') = cos{u'^) and 
a^{v}.,v?) = v}sin{u^), we can compute 



, i , I 2\\ ( —sin(u^) 

which satisfies {H2) (under assumptions (i) and (ii)). Therefor Theorem \l.l\ gives the 
existence of a solution to lil.6\} with (i) and (ii). 

Based on the same type of gradient entropy inequality p.3p , it was proved in Cannone 
et al. [Sj the existence of a solution in the distributional sense for a two-dimensional 
system of two transport equations, where the velocity vector field is non-local. 

The uniqueness of the solution is strongly related to the existence of regular (Lipschitz) 
solutions (see Theorem I7.7p . Let us remark that equation (lPll- (IIDI ) does not create 
shocks because the solution (given in Theorem ll.ip is continuous. In this situation, it 
seems very natural to expect the uniqueness of the solution. Indeed the notion of en- 
tropy solution (in particular designed to deal with the discontinuities of weak solutions) 
does not seem so helpful in this context. Nevertheless the uniqueness of the solution is 
an open problem in general (even for such a simple system). 

We ask the following Open question: 
Is there uniqueness of the solution given in Theorem 11.11 ? 

Now we give the following existence and uniqueness result in [14^^'°°([0, T) x M)]^^, in a 
special case to simplify the presentation. More precisely we assume 

{HV) a\u) = Y^ A,jV? for i = 1, . . . , M and for all u^U, 

i=l,...,M 



(if 2') Y, ^ni^i, > for every ^ = (^i, -,6/) e [0, +oo; 

i,j=l,...,M 



M 



Theorem 1.5 (Existence and uniqueness of W^'°° solution for a particular 

Assume {HI'). For T > and all nondecreasing initial data Uq G [iy^'°°(M)]*^, the sys- 
tem ^-^^) admits a unique solution u G [W^^'°°([0, T) x M)] , in the following cases: 

i) M>2 and Aij > 0, for all j > i. 

ii) M >2 and Aij < 0, for all i -/- j and {H2'). And then for all (t,x) G [0,T) x M wje 

have 

Y 9.u\t,x) < sup J2 9.ul{y). (1.7) 



--1,...,M ''^"" i=l,...,M 



Remark 1.6 (Case of M = 2) 

In particular for M = 2, if {HI'), {H2') and {H3) satisfied then we have, by Theorem 
\1.5\ the existence and uniqueness of a solution in [W^'°°{[0, T) x M)] of ^)-^^^. 

In these particular cases of the matrix A, we can prove that d^u^ for i = 1, . . . , M, are 
bounded on [0, T) x M. Thanks to this better estimates on d^u^, and then on the velocity 
vector field Au^ we prove here the uniqueness of the solution. 

In the case of the matrix ^ = I ^ ^ J , it was proved in El Hajj, Forcadel [E], the 



existence and uniqueness of a Lipschitz viscosity solution, and in A. El Hajj [15], the 
existence and uniqueness of a strong solution in W^^^ ([0,T) x R). 



1.3 Application to diagonalizable systems 

Let us first consider a smooth function u = {u^, . . . ,m*^), solution of the following non- 
conservative hyperbolic system: 

dtu{t, x) + F{u)d^u{t, x) = 0, u{t, x) eU, xeR, te (0, T), 

(1.8) 
u{x, 0) = uo{x) X eM, 

where the space of states U is now an open subset of M*^, and for each m, F{u) is a 
M X M-matrix and the map F is of class C^{U). We assume that F{u) has M real 
eigenvalues Ai(m), . . . , Am(^), and we suppose that we can select bases of right and left 
eigenvectors ri{u), Uiu) normalized so that 

|rj| = 1 and /j ■ r^ = 5ij 



Remark 1.7 (Riemann invariant) 

Recall that locally a necessary and sufficient condition to write 

li{u) = Vu^iiu), 

is the Frobenius condition li A dli = . In that case the function (pi{u) is solution of the 
following equation 

{iPi{u))t + \i{u){ipi{u))x = 0. 
We recall that then (pi{u) is called a i- Riemann invariant (see Sevennec J3^ and Serre 



fSdj Vol II])). If this is true for any i, we say that the system lll.8\) is diagonalizable. 
Our theory is naturally applicable to this more general class of systems. 
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1.4 A brief review of some related literature 

Now we recall some well known results for system (11 .Sp . 

For a scalar conservation law, this corresponds in (11 .Sp to the case M = 1 and 
F{u) = h'{u) is the derivative of some flux function h, the global existence and unique- 
ness of BV solution established by Oleinik [34j in one space dimension. The famous 
paper of Kruzhkov [28] covers the more general class of L°° solutions, in several space 
dimension. For another alternative approach based on the notion of entropy process 
solutions, see Eymard et al. [1^, see also the kinetic formulation P. L. Lions et al. |33j . 

We now recall some well-known results for a class of 2 x 2 strictly hyperbolic systems n 
one space dimension. Here i.e F{u) has two real, distinct eigenvalues 

Ai(m) < \2{u). 

Lax [30] proved the existence and uniqueness of nondecreasing and smooth solutions of 
the 2x2 strictly hyperbolic systems. Also in case of 2 x 2 strictly hyperbolic systems 
DiPerna [12], [13] showed the global existence of a L°° solution. The proof of DiPerna 
relies on a compensated compactness argument, based on the representation of the weak 
limit in terms of Young measures, which must reduce to a Dirac mass due to the presence 
of a large family of entropies. This results is based on the idea of Tartar [39] . 



For general M x M strictly hyperbolic systems; i. e. where F{u) has M real, distinct 
eigenvalues 

Ai(m) <■■■< Aa/(m), (1.9) 

Bianchini and Bressan proved in [6] a striking global existence and uniqueness result of 
BV solutions to system (11.81) . assuming that the initial data has small total variation. 
Their existence result is a generalization of Glimm result [20], proved in the conserva- 
tion case; i.e. F{u) = Dh{u) is the Jacobin of some flux function h and generalized by 
LeFloch and Liu [31], [32] in the non-conservative case. 



We can also mention that, our system (|P|) is related to other similar models, such as 
scalar transport equations based on vector flelds with low regularity. Such equations 
were for instance studied by Diperna and Lions in [H]. They have proved the existence 
(and uniqueness) of a solution (in the renormalized sense), for given vector fields in 
L^((0,+oo);iy/;J(M^)) whose divergence is in L^((0, +oo); L°°(M^)). This study was 
generalized by Ambrosio [2], who considered vector fields in L^ ((0, +oo); i?V/oc(K^)) 
with bounded divergence. In the present paper, we work in dimension A^ = 1 and prove 
the existence (and some uniqueness results) of solutions of the system (|P]l-( lIDp with a 
velocity vector field a^{u), z = 1, . . . , M. Here, in Theorem 11.11 the divergence of our 
vector field is only in L°°((0, +oo), LlogL(M)). In this case we proved the existence 
result thanks to the gradient entropy estimate (11.31) . which gives a better estimate on 



the solution. However, in Theorem 11.51 the divergence of our vector field is bounded, 
which allows us to get a uniqueness result for the non-linear system f |Pl) . 



We also refer to Ishii, Koike [25] and Ishii [24], who showed existence and uniqueness of 
viscosity solutions for Hamilton-Jacobi systems of the form: 

dtu' + Hi{u, Du') = with u = {u% € M^, for xeR^, te (0, T), 

(1.10) 

u^{x, 0) = ul{x) X eR, 

where the Hamiltonian Hi is quasi-monotone in u (see Ishii, Koike [25], Th.4.7]). This 
does not cover our study since our Hamiltonian is not necessarily quasi-monotone. 



For hyperbolic and symmetric systems, Garding has proved in [18] a local existence and 
uniqueness result in C{[0,T);H'(R^))nC^{[0,Ty, H'-^R^)), with s > f + 1 (see also 
Serre [36l Vol I, Th 3.6.1]), this result being only local in time, even in dimension N = 1. 

1.5 ]V[iscellaneous extensions to explore in a futur work 

1. In Theorem [HT] we have considered the study of a particular system only to simplify 
the presentation. This result could be generalized to the following system 

dtu' + a\u,x,t)d^u' = h\u,x,t) on (0, T) x M and for t = l,...,M, (P') 

with suitable conditions on a* and h\ 



2. If we consider the case where the system (|P|) is strictly hyperbolic. Based in the 
result of Bianchini, Bressan [6], we could also prove the uniqueness of the solution, 
whose existence is given by Theorem 11.11 

3. We could also extend Theorem 11.51 to system (iPj) . where we replace (i) and (ii) by 
the following condition 

i') For M > 2, a^iu, x,t)>0 for j > i and for all {u,x,t) eU xRx [0, T). 
ii') For M > 2, 

a'j{u, x,t) <0 for all {u,x,t) e U xR x [0, +oo), for all i ^ j, 
and we assume that for any Vi G R^,Xi G R, the matrix 

bij{t) = a'j{vi,Xi,t) 
satisfies for allt > 

iH2") Yl %We4>0 for all ^ = (6, -.eM) G [0, +oo)^. 

i,j=l,...,M 

8 



4. We could also prove the uniqueness result in case of W^'°° solution among weak 
solution, (and in particular any weak solution is a viscosity solution in the sense of 
Crandall-Lions [inilll]). 

5. We could propose a numerical scheme and try to prove its convergence. 

6. Applications to other equations: Euler, p-systems. 



1.6 Organization of the paper 

This paper is organized as follows: in the Section 2, we approximate the system (lP|) and 
the initial conditions. Then we prove a local in time existence for this approximated 
system. In Section 3, we prove the global in time existence for the approximated system. 
In the Section 4, we prove that the obtained solutions are regular and non-decreasing 
in X for all t G (0,T). In Section 5, we prove the gradient entropy inequality and some 
other e-uniform a priori estimates. In Section 6, we prove the main result (Theorem ll.il) 
passing to the limit as e goes to and using some compactness properties inherited from 
our entropy gradient inequality and the a priori estimates. In Section 7 we prove some 
uniqueness results in particular cases (Theorem II. 5p . An application to the dynamics of 
dislocation densities given in Section 8. Finally, in the Appendix, we recall the proof of 
uniqueness of Lipschitz solution to system ([P]). 



2 Local existence of an approximated system 



The system (|P|) can be written as: 

dfU + a{u) o dxU = 0, 



(2.11) 



where u := («*)!,..., a/, «'(«) = (ct*('ii))i,...,M and UoV is the "component by component 
product" of the two vectors U,V E M^. This is the vector in M^^ whose coordinates are 
gWenhyiUoV)r-=U,V,: 



U2 





V2 
_Vm _ 


= 


U2V2 

_ UmVm _ 



Now, we consider the system (12.111) . modified by the term ed^xU, where d^ 



dx- 



■, and 



for smoothed data. This modification brings us to study, for all < £ < 1, the following 
system: 



dtu^ — edxxu'^ = —a{u'^) o dxu'^ ■, 



(Pe 



with the smooth initial data: 

ti^(x, 0) = ul{x), with Uq{x) := Uq * "q^ix), (-^-^e) 

where % is a moUifier verify, %{■) = -??(-), such that r] G C^(]R) is a non-negative 
function and J^t] = 1. 

Remark 2.1 

By classical properties of the moUifier {f]e)e o.nd the fact that Uq G [L°°(M)]*^, then 
ul G [C°^(M)]*^ n [W"''°°{R)]^ for all meN. 



The global existence of smooth solution of the system (j^j) is standard. Here, we prove 
this results only to ensure the reader. 

The following theorem is a local existence result (in the "Mild" sense) of the regularized 



system (|^-( |JD^ . This result is achieved in a super-critical space. Here particularly 



we chose the space of functions [C([0,T);X(]R))] , where 

X{R) = {ue L°°{R) such that d^u G L^{R)}. (2.12) 

This space is a Banach space supplemented with the following norm 

||'"IU(M) = ||'"||l°°(R) + ||f^x^||L8(R)- 

Here the espace Lp(M) with p = 8 will simplify later in Lemma [4T] the Bootstrap argu- 
ment to get smooth solution. 

In this Section, we will prove the following 

Theorem 2.2 (Local existence result) 

For all initial data u^ G [X(]R)]^^ there exists 

T* = T*{Mo,e) >0, 
such that the system [P^- lJIDeD admits solutions m^ G [C([0, T*);X(R))] . 



In order to do the proof of Theorem 12.21 in Subsection 12.21 we need to recall in the 
following Subsection some known results. 

2.1 Useful results 

Lemma 2.3 (Mild solution) 

Let T > 0, and u^ G 
with u^{t) = u'^{t, ■): 



Let T > 0, and u^ G [C([0,T); X(]R))] he a solution of the following integral problem 



U'it) = Se{t)ul - I Seit - S) HU^S)) O d^,u' (s)) ds, {I N^ 

Jo 
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where Ss{t) = Si{et) such that Si{t) = e*^ is the heat semi-group. Then m^ is a solution 
of the system [P^- fjIDsD in the sense of distributions. 



For the proof of this lemma, we refer to Pazy [35l Th 5.2. Page 146]. 

Lemma 2.4 (Picard Fixed Point Theorem, see 126^) 

Let E be a Banach space, let B : E x E — > E be a continuous map such that: 

\\B{x,y)\\E <r]\\y\\E for all x,y e E, 
where r] is a positive given constant. Then, for every Xq G E, if 

<r] <1, 

the equation x = Xq + B{x,x) admits a solution in E. 



In order to show the local existence of a solution for ( [JA^eD , we will apply Lemma 12.41 in 
the space E = [L°°((0, T); X(M))]*^ 

Lemma 2.5 (Time continuity) 

LetT > 0. If u^ e [L°°{{0,T);W^'P{R))]^^, 1 < p < +oo, are solutions of integral 
problem ^TTQ , then u' G [C{[0,T);W^'P{R))]^ . 

For the proof of Lemma [231 see A. Pazy [35l 7.3, Page 212]. 

Lemma 2.6 (Semi-group estimates) 

Let I < p < q < +00. Then for all f G Lp(M) and for all t > 0, we have the following 
estimates: 






LP( 



^) \\Ss{t)f\\LHm<ct--'-^-^ 

where C = C{e) is a positive constant depending on e. 

For the proof of this Lemma, see Pazy [35l Lemma 1.1.8, Th 6.4.5]. 

2.2 Proof of Theorem [2721 



Our goal is to show local existence of a solution of IjP^ using the Picard fixed point 
Theorem. To be done according Lemma [2731 it is enough to prove the local existence for 
the following equation: 

u'{t) = S,{t)ul - [ S,{t - s) {a{u%s)) o d,u%s)) ds, 

■^0 (2.13) 

= SMnl + B{u',u%t), 
11 



with B{u,v){t) = — Sei^t — s) {a{u){s) odxv{s))ds. 
Jo 

If we estimate B{u,v), we will obtain, for all u,v & [L°°((0,T);X(M))]*^, where X{ 
defined in (12.121) . the following: 



\B{u,v)m 



[nmv 



+ 



Se{t - s) {a{u{s)) o dxv{s)) ds, 



dr,Se{t - s) {a{u{s)) o dj:v{s)) ds, 



[L»{R)]^ 



where for a function / = (/^, . . . , f^^) G [X(R)]^''^, we note here 



[x(M)]A/ = sup ||/1l^(m)+ sup \\dxf 

i=l,...,M i=l,...,M 



(2.14) 



Using Lemma [2^61 (i) with p = 8, g = oo for the first term and Lemma [221 (ii) with p = 8 
for the second term, we obtain that : 



B{u,v){t)\\ix{R)]M< C -^||a(M(s))9^f(s)||[^2(^)jMrfs, 



(t-s) 



* 1 



-\\a{u{s))dxv{s)l^s(mMds 



/o (t-s) 
We use the Holder inequality, and get, for all < T < 1: 

< CT^ \\v 



(2.15) 



L°°((0,T);X(R))]" ) 

where C{MQ,e). Moreover, we know by classical properties of heat semi-group (see A. 
Pazy [35]): 



\Ssi't)ul\\[L°°{(0,T);X{R))]M < \\ul\\[x(R)]M ■ 



Now, taking 



we can easily verify that 



(ni=min(^,l 



C(T*)2 < 1. 



(2.16) 



(2.17) 
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By applying the Picard Fixed Point Theorem (LemmaEH) with E = [L°°((0, T*); X(M))]^, 
this proves the existence of a solution u^ E [L°°((0,T*); X(M))]*^ for ^A . 



Then, according to Lemma [275l we deduce that the solution is indeed in [C([0, T*); X 
This proves, by Lemma [231 the existence of a solution in [C([0, T*); X(R))]*^, which 



M 



satisfies the system ([^-( |/iJ^ in the sense of distributions. □ 



3 Global existence of the solutions of the approxi- 
mated system 



In this Section, we will prove the global existence of solution for the system ([^])-( |/iJ^ . 
Before going into the proof, we need the following lemma. 

Lemma 3.1 {L°° bound) 

Let T > 0. Ifu^e [C([0,T);X(M))]^^ ts a solution of system (F^-^D^) with initial 
data Uq G X(M), then 

11"" ll[L°°([0,T)xM)]*^ — IPoll[L°o(IR)]*^ 

The proof of this Lemma is a direct application of the Maximum Principle Theorem for 
parabolic equations (see Gilbarg-Trudinger [19l Th.3.1]). 

Remark 3.2 

Thanks to the previous Lemma, we notice that we can take the box U defined in lll.l\) as 
the following 

TT TT-A/ r II £-*ll II £1*11 1 

U =iij=i[-||Mo I|l°°(ir), 11^0 IU°°(i:)J- 
For fixed e, this definition guarantee that Mq do not change in the course of time. 



The result of this Section is the following. 

Theorem 3.3 (Global existence) 

Let T > and < £ < 1. For initial data u^ E [X{M)] satisfying (HI) and {H2). 
Then the system (F^-(TD^, admits a solution u^ e [C([0,T);X(M))]^^ with u^{t,-) 
satisfying {HI) and {H2) for all t e (0, T). Moreover, for all t e (0,T), we have the 
following inequalities: 

lk^''(^>-)IU°°(iR) < lko'*IU°=(K)' fori = l,...,M, (3.18) 

Proof of Theorem ISTSt 

We are going to prove that local in time solutions obtained by Theorem 12.21 can be 
extended to global solutions for the same system. 
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We argue by contradiction: assume that there exists a maximum time T^ax such that, 



we have the existence of solutions of the system ([^])-( |/Z)^ in the function space 

[C([0,T^„,);X(R))]*^. 



For every small enough 6 > 0, we consider the system ([^ with the initial condition 



From Theorem 12.21 to deduce that there exists a time T*(Mo,e), independent of S (see 
Remark [3^ . such that the system d^j ) with initial data Uq has a solution u'^'^ on the 
time interval [0,T*). Then for 

^0 = [Tmax — S) +T , 

we extend u'^ on the time interval [0, Tq) as follows, 



u'^{t,x) 



U'{t,x), for t e [0,Tmax - S], 



M^' (t,x), for t G [Tmax " S,Tf^ 



0) 



and we can check that u'^ is a solution of ([^|)-( |/iJ^ on the time interval [0,Tq). But 
from Lemma (13 .ip we know that the time T* is independent of 5 (see Remark l3.2l) . which 
implies that Tq > Tmax and so a contradiction. 
The inequalities (I3.18P is a consequence of Lemma [3Tl □ 

4 Properties of the solutions of the approximated sys- 
tem 



In this section, we are going to prove that the solution of ([^1)-( |/iJ^ obtained by Theorem 
12.21 is smooth and monotone. 

Lemma 4.1 (Smoothness of the solution) 

Let T > 0. For all initial data m§ G [X(M)]^^ where d^uf^ G [^^'"'^(M)]^^ for all meN, 
1 < p < +00. 

If u^ is a solution of the system JP^- IJID^ ), such that u^ G [C([0,T);X(M))] and 
d^u' G [L°°((0,T);Li(R))]*^ then u' G [C°°([0,T) x R)]^^ and satisfies, 

u' e [W"''P{{0, T) X R)]^^ , for all 1< p < +oo and m eN\ {0}, (4.19) 

Proof of Lemma 14.11 

Step 1 (Initialization of the Bootstrap): 

For the sake of simplicity, we will set 
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F[u^] = —a{u^) o dxU^. 

From the fact that u^ G [C([0,T);X(M))]^^ and d^u^ e [L°°((0,T); L\R))f^, we deduce 
that d^u', F[u'] e [L^{iO,T) x R)]^ n[L^{{0,T) x M)]^, which proves by interpolation 
that 

d^u', F[u'] e [LP{{0,T) X R)]^ for all 1 < p < 8. (4.20) 

Because u'^ is a solution of ( [^1 ), we see that 



dtu' - ed,,u' = F[u% (4.21) 

a^^u^ - 6d,,,u' = d,F[u% (4.22) 

Applaying the classical regularity theory of heat equations on (14.211) . we deduce that: 

dtu' and d^^u' G [LP{{0,T) x M)]^^ , for all 1< p < 8. (4.23) 

For more details, see Ladyzenskaja [2Sl Theorem 9.1]. But we know that 

d^F[u'] = -a{u') o d.^^u' - Da{u')d^u' o d^^u' (4.24) 

We notice that thanks to this better regularity on u'^ ( (14.201) and (14.231) . and by the 
Holder inequality we can easily prove that 

da:F[u'] G [LP((0, T) X R)f^ for all 1< p < 4. 

Now, we apply again the classical regularity theory of heat equations on (|4.22l) . to deduce 
that: 

dt^u' and d^^^u' G [LP((0,T) x R)f^ , for all 1< p < 4. (4.25) 

We know that 

dtF[u'] = -a{u') o dt^u' - Da{u')dtu' o d^u' (4.26) 

Thanks this previous regularity on u^ , we obtain by the Holder inequality that 

dtF[u'] G [LP((0, T) X R)]*^ for all 1< p < 4. 
Which gives that 

a^u^ F[u'] G [iy^'P((0,T) X R)]*^ for all 1< p < 4, 

and by the Sobolev embedding that d^u' G [Lp((0,T) x R)]*^ for all 1 < p < oo. 
Step 2 (Recurrence): 
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(Hm) 



Now, we use the same steps, we can prove by recurrence that for all m G N if, 

d^u', F[u'] G [IV^'PiiO, T) X R)]^ for all 1 < p < 4, 
then 

Indeed, as in (14.231) we can deduce here that 



M 



dtu' and d^^u' G [Vr"''P((0, T) x M)Y' , for all 1< j9 < 4 



(4.27) 



and From (14.241) . because 9^m^ G [L°°((0,T) x M)] , we can obtain here that 

d^F[u'] G [Vr'"'P((0, T) X M)]^ for all 1< p < 4. 
Which proves that, as in (14.251) that 

dt^u' and d^^^u' G [iy"'P((0, T) x M)]^ , for all 1< p < 4, 
and From (14.261) . we deduce that 



(4.28) 



dtF[u^] G [iy'"'P((0,T) X 



for all 1 < p < 4, 



and then 



M 



d^u', F[u'] G [W^"^+^'P((0,T) X R)J^" for all 1< p < 4, 
Which proves by the Sobolev embedding the results. 



D 



Lemma 4.2 (Classical Maximum Principle) 

Let T > 0. For all initial data m§ G [X(M)]*^, where 5^.Mg G [Vr™'P(M)]^ /or a// m G N, 
1 < p < +CXD, and satisfying {H3). 



If u'^ is a solution of the system lP^- lj^IDir\) , such that u^ G [C([0,T);X 



M 



and 



d^u' G [L°°((0, T); Li(R))]^, t/ien we have for i = 1, . . . ,M , d^u''' > on (0, T) x 
Proof of Lemma 14.21 



We first derive with respect to x the system ([^-( |/iJ^ , and get for w"^ = (w^'*)j=i,..., 



,,M 



dxU'^ 



dtw"^ - edxxw'^ + a{u'^) o d^w'^ + Da{u)w'^ ow'^ = 0. 



tM 



Since u"^ G [C°°([0, T) x M)] , we see, for i = 1, . . . , M, that w^'^ is smooth and satisfies 
w^'\0,x) = dxUQ^ > 0. From the classical maximum principle we deduce that w^'"^ > 
on[0,T)xM. D 
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Remark 4.3 (L^ uniform estimate on 9xM^) 

Because d^u^'"^ > 0, for i = 1, . . . , M , we deduce from Lemma lSAl that: 

\\dxU ||[Loo([o,T);Li(IK))]" — ^11"" \\lL°°{[0,T)xR)]'^' — '^W'^oW [L°° (R)]'^' ■ (4.29) 



j=l,...,M 



Corollary 4.4 (global existence of nondecreasing smooth solutions) 

Let T > 0. The solution given in Theorem \2.S\ can be chosen such that u'^ = {u 
smooth, satisfies ^4-19^ and for each 2 = 1,..., M , d^u'^'^ > on (0, T) x M. 

The proof of Corollary 14.41 is a consequence of Theorem 12.21 and Lemmata 14.11 14.21 and 
Remark 14.31 

5 e-Uniform a priori estimates 



In this Section, we show some 5-uniform estimates on the solutions of the system l]P^ 



d/ZJgl) . These estimates will be used in Section [6] for the passage to the limit as e tends 



to zero. 



Lemma 5.1 {L°° bound on u^ and L^ bound on S^m^) 



M 



Let T>0, 0<e<l and function Uq G [L°°(R)] satisfying {H3). Then the solution 



of the system [F^ - ^TD^ ) given in Theorem \3. Si with initial data Uq = Uq * rj,;, satisfies 



the following e-uniform estimates: 

{El) ||lt^||[Loo((o,T)xR)]'*^ — ll^o||[L°o(R)]*^ ! 

\E2) ||c^a;«^||[ioo((o_T),Ll(IR))]^-' — 2 |Po||[Lo>=(K)]M , 

Proof of Lemma 15.1b 

First, we remark that lid^UQ > 0, then d^u^ = {O^Uq) *rj^{x) > (because rj is positive). 

The fact that uq e [L°°(M)]^ and d^uo > 0, we obtain that d^uo G [L\R)f^. 

By classical properties of the mollifier (r/e)^ we know that if Uq G [L°°(M)] and 

d^uo G [L\R)]^^ we have m§ G [X(M)]^^ and d^u^ G [W'^'PiR)]^ for all m G N, 

1 < P < +00. 

Now, we use Lemma 13.11 and Remark 14.31 we deduce by the classical properties of the 
mollifier {El) and {E2). 

Before going into the proof of the gradient entropy inequality defined in f l5.30p . we 
announce the main idea of this new gradient entropy estimate. Now, let us set for w > 
the entropy function 

f{w) = wlnw. 
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For any non-negative test function ip G C^(IR x [0,+oo)), let us define the following 
"gradient entropy" with w* := d^u^: 

■^^ V=1,...,M / 

It is very natural to introduce such quantity N{t) which in the case </? = !, appears to be 
nothing else than the total entropy of the system of M type of particles of non-negative 
densities w\ Then it is formally possible to deduce from ([P]) the equality in the following 
new gradient entropy inequality for all t > 

^(^) + [ v( y2 «>'^M d^ < Rit) for t > 0, (5.30) 

"^^ ^^ Vj=i,...,M ' / 

with the rest 



-^^ y \j=i,...,M / \i=i,...,M 



w*) > dx. 



where we only show the dependence on t in the integrals. We remark in particular that 
this rest is formally equal to zero if ip = 1. 

To guarantee the existence of continuous solutions, we assumed in {H2) a sign on the 
left hand side of inequality (|5.30l) . 

For we return this previous calculate more rigorous, we prove actually the following 
gradient entropy inequality 

Proposition 5.2 (Gradient entropy inequality) 

Let T>0, 0<£:<1 and function uq G [L°°(]R)] satisfying {H3). We consider the 
solution u^ of the system ^P^- ^TD^ given in Theorem \3.3\ with initial data Uq = Uo*ri^,. 
Then, there exists a constant C{T,M,Mi, ||-Uo||[loo(r)]A'/, ||t^x'Wo||[LiogL(R)]A-^ such that 

^^*)+/ / E aij{u')w'''w''^ <C, with N{t)= f J2 f{w''')dx. (5.31) 

where w^ = (if^'*)j=i,...,Af = dxU'^ and f is defined in (jJ.-^l ). 
For the proof of Proposition 15.21 we need the following Lemma: 

Lemma 5.3 (LlogL Estimate) 

Let {ris)s be a non-negative mollifier, f is the function defined in (jJ.-^D and h G L^(M) is 
a non-negative function. Then 
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i) / f{h) < +00 if and only if h G LlogL(R). 

Jr 

a) If h E LlogL(R) the function h^ = h * i]^ E LlogL{R) satisfies 

\\h- hs\\LiogL(R) ^ as e^O. 

The proof of (i) is trivial, for the proof of (ii) see R. A. Adams [H Th 8.20] for the proof 

of this Lemma. 

Proof of Proposition 15.2b 

Remark first that the quantity N{t) is well-defined because w'^ E [L°°((0, T); L^(R))] n 
[L°°((0,T); L^(]R))] (by Theorem 12.21 and Corollary 14. 4p and we have the general in- 
equality ^ < wlogw < w"^ for all w > 0. 

From Theorem 14.41 we know that w^'* and smooth non-negative function. Now, we 
derive N{t) with respect to t, this is well-defined because for i = 1, . . . , M, we have 

< e||M;"'^||ioo((o,T);Li(R)) and for all m G N, w'-^' E PF"'°°((0,T) x M) (see fl419D V 

— e 

Finally, we get that. 



l-« 



-'^1=1,.. .,M 
•^^i=l,...,M 



Ji J2 



-^ V . ^ 






w'^') 



But, it is easy to check that 

/■// \ _ / lii(^) + 1 if X > 1/e, , „„, \ _ / ^ if X > 1/e, 

^ ^^^ ~ \ if < X < 1/e, and / (xj - j ^ if < x < 1/e. 

This proves that J2 < 0. To control Ji, we rewrite it under the following form 

Ji= I Y ci\u')9\w'^')d,w^-\ 



where 



_ ^ ~ \ if X > 1/e, 
^^^^ ~ "i ' if < X < 1/e, 
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Then, we deduce that 



■^1 = / E ^\u')U9{w'')) 

■^^ i=l,...,M 

'^ M ■ ■_-| nj 



i,j=l,...,M 

Jll Jl2 
^N . . ^^ 



'^''^ i,j=l,...,Af "^'^i,j=l,...,M 

From (-ff2), we know that Jn < 0. We use the fact that \g{x) —x\<\ for all x > and 
(ifl), to deduce that 

\Jl2\ < -M Ml |k^'1|[j:oc((o,T),Li(K))]*^ 

< -M Mi||no||[L«>(R)]M 

where we have use Lemma EH] (-E2) in the last line. Finally, we deduce that, there exists 
a positive constant C{\\uq\\[ioc^-^-j^m , Mi, M) independent of e such that 

-N{t) <Jll + Jl2 + J2 

at 

< Jll + c. 



Integrating in time we get by Lemma 15.31 there exists a another positive constant 

C(T, M, Ml, ||mo||[l-'(r)]m, \\dxUo\\[L\ogL{m.)]M) independent of e such that 



+ /7 S <'< 

Jo J^i^j=l^,,,M 



^(t) +11 > ^ a[^{u')w''^w''' <CT + N{0) < C. 

D 

Lemma 5.4 (W^^'^ estimate on the time derivatives of the solutions) 

Let T>0,0<e<l and function uq G [L°°(R)] satisfying {H3). Then the solution 
of the system ^P^ - [TD2 ) given in Theorem \3. Si with initial data u^ = Uq * rj^, satisfies 
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the following e-uniform estimates: 

where W~^^^{R) is the dual of the space W^'^{R). 

Proof of Lemma 15.4b 

The idea to bound dtu'^ is simply to use the available bounds on the right hand side of 



the equation [P^] 
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We will give a proof by duality. We multiply the equation ([^ by G [^^((O, T), iyi'°°(]R))] 
and integrate on (0, T) x M, which gives 



M 



f 

J(0. 



h 



h 






(f) dtu^ = ^ / 't> d^x^"^ ~ 4> a{u^) o dxU^ ■ 

'(0,T)x]R i(0,r)xR i(0,T)xR 

We integrate by parts the term Ji, and obtain that for < e < 1: 



\h\< 



dx4>dxu'^ 



(0,T)xI 



< 7'||aa;0||j^2((o_T),L°°(M))]*Ml^a;M II [l2((0,T),L1(R))]" ' 



< 2T2||0||j^2((o,t),m/i.°°(r))]*mI'"o||[lo>=(ir)]M, 
here, we have used the inequality 



\\dxu' 



< 2r2||Mo| 



[2,oo(]g)]M, 



(5.32) 



(5.33) 



l[L2([0,T);Li(M))]^ 

which follows from estimate (14.290 for bounded and nondecreasing function m^. Similarly, 
for the term I2, we have: 



I/2I <Mo||m| 



[L°°((0,T)xI 



[L2((o,T),L° 



l<9:,:M^ 



l[L2((o,T),L3 



im^M*^; 



< 2T^Mo||mo|'^ 



(5.34) 



[L°°(R)]*^ll'^ll[L2({0,T),iyi. 



Finally, collecting (15.321) and (15.34p . we get that there exists a constant C = C(T, Mq) 
independent of < £ < 1 such that: 



(pdtu' 



(0,T)xI 



<C 1 



\uo\ 



[L° 



[L2({o,T),VKi 



which gives the announced result where we use that L {{0, T), W ' (M)) is the dual of 
L2((0,T), Vr^'°°(R)) (see Cazenave and Haraux P Th 1.4.19, Page 17]). □ 

Corollary 5.5 (e-Uniform estimates) 

LetT > 0,0 <e <1 and function Uq G [L°°(M)]*^ satisfying (HI) and {H2). Then the 
solution of the system /f^ -i f/DTP given in Theorem \3.S\ with initial data ul = uq * %, 
satisfies the following e-uniform estimates: 



\\OxU ||[Loo((o_T);LlogL 

where C = C{T, M, Mq, Mi \\uol^ 



*^ + 11"" ll[L°=((0,T)xIR)]*^ 



M , ||<9xMo||rLiogL(M)lAO- 



<c. 



We can easily prove this Corollary collecting Lemmata 15. 11 15.4] and 15.31 and Proposition 

El 
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6 Passage to the limit and the proof of Theorem 11.1 



In this section, we prove that the system (lPll- (IIDp admits solutions u in the distributional 
sense. They are the limits of m^ given by Theorem 13.31 when e — > 0. To do this, we will 
justify the passage to the limit as e tends to in the system ([^])-( |/JJ^ by using some 
compactness tools that are presented in a first Subsection. 

6.1 Preliminary results 

First, for all / open interval of M, we denote by 

L\ogL{I) == lfeL\l) such that /" |/| In (1 + |/|) < +oo 

Lemma 6.1 (Compact embedding) 

Let I an open and bounded interval of M. // we denote by 

y^i,LiogL^j) = {m G L\I) such that d^u G LlogL(/)}. 
Then the following injection: 

is compact. 

For the proof of this Lemma see R. A. Adams [H Th 8.32]. 

Lemma 6.2 (Simon's Lemma) 

Let X , B , Y be three Banach spaces, such that 

X ^-^ B with compact embedding and B "—^ Y with continuous embedding. 
Let T > 0. If ('U^)^ is a sequence such that, 

||'"^||l°°((0,T);X) + ||'"^||l°°((0,T);B) + l|f^t^^ llL9((0,T);y) - C*, 

where q > 1 and C is a constant independent of e, then ('U^)e is relatively compact in 
C((0,T);i?). 

For the proof, see J. Simon [38], Corollary 4, Page 85]. 

In order to show the existence of solution system (|Pj in Subsection [621 we will apply this 
lemma to each scalar component in the particular case where X = W^'^°^{I), B = L°°{I) 
and Y = W-^'\I) := {W^'^{I)y. 

We denote by Kexp{I) the class of all measurable function m, defined on J, for which, 
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/"(el^l-l) <+oo. 



The space EXP{I) is defined to be the linear hull oi Kexp{I)- This space is supplemented 
with the following Luxemburg norm (see Adams [1] (13), Page 234] ): 



\u\\EXPii) = inf <^ A > : / (^e'5^' - l) < 1 I 



Let us recall some useful properties of this space. 

Lemma 6.3 (Weak star topology in LlogL) 

Let Eexp{I) be the closure in EXP{I) of the space of functions bounded on /. Then 

Eexp{I) is a separable Banach space which verifies, 

i) LlogL(J) is the dual space of Eexp{I)- 

a) L^{i)^E,xpii). 

For the proof, see Adams P Th 8.16, 8.18, 8.20]. 

Lem,m,a 6.4 (Generalized Holder inequality, Adams [H, 8.11, Page 234]) 

Let f E EXP{I) and g G LlogL(J). Then fg E L\l), with 

||/fi'||Li(/) < 2||/||£;xP(/)||fi'||LlogL(/)- 

The following Lemma, we allow to define later the restriction of a function / E W~^'^(R) 
on all open interval / of M. 

Lemm,a 6.5 (Extension) 

For all open interval I of M, there exists a linear and continuous operator of extension 
P : W^'°°{I) -^ l^i'~(M) such that 

i) Pu\^ =uforuE iyi'°°(/). 

ii) ||PM||vyi,oo(K) < ||m||p^/i,oo(/) foru E iy^'°°(J). 

for the proof of this Lemma see for instance Brezis [7l Th.8.5]. 



Thanks this Lemma, we can notice that, \i f E W ^'H^^)? where W ^'^( 
(iy^'°°(M))', we can define, for all open interval I of M, the function f\j as the following 

< f\iih >w-^'^{i),w^'°°{i)=< fjPh >vK-i^i(K):VKi.°°(M) • 
for all /le W^^°^{I). 
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6.2 Proof of Theorem [TTI 
Step 1 (Existence): 

First, by Corollary 15.51 we know that for any T > 0, the solutions u'^ of the sys- 
tem ([^-( |/D^ obtained with the help of Theorem 13.31 are e-uniformly bounded in 
[L°°((0, T) X M)] . Hence, as e goes to zero, we can extract a subsequence still denoted 
by M^, that converges weakly-* in [L°°((0,T) x R)] to some limit u. Then we want to 
show that u is a solution of the system (lPl)- (IIDI ). Indeed, since the passage to the limit 
in the linear terms is trivial in [P'((0,T) x R)] , it suffices to pass to the limit in the 
non-linear term. 

According to Corollary 15.51 we know that for all open and bounded interval / of R there 
exists a constant C independent on e such that: 



ll"'ll[L-{(0,T);H/i.^l°g^(/))]" + 11^11 [L-((0,T)x/)]*^ + ll^««^ll [L2((0,T);H/-i.1(/))]« < C' 

From the compactness of iy^'^'°§^(J) ^^ L'^(T) (see Lemma 16731 (i)), we can apply 
Simon's Lemma (i.e. Lemma Q, with X = [W^'^^"^^{I)Y\ B = [L°^(/)]^^ and 
Y = [W-^'^{I)]^, which shows that 



u' is relatively compact in in [L°°((0,T) x I)f^ ^ [L\{0,Ty, L'^{I))]^ . (6.35) 

Then form continuous injection of L°°{I) --^ Eexp{I) (see Lemma [6731 (ii)), we deduce 
that, 

u^ is relatively compact in [L^{{0,T);Eexp{^))] ■ (6.36) 

On the other hand, by Corollary 15.51 we notice that dxU'^ is e-uniformly bounded in 
[L~((0,r);LlogL(/))]^^. Moreover, the space [L°°((0,T); LlogL(/))]^^ is the dual 
space of [L^{{0,T);Eexp{I))] , because LlogL(/) is the dual space of Eexp{I) (see 
Lemma 16731 (ii) and Cazenave, Haraux [91 Th 1.4.19, Page 17]). Then, up to a sub- 
sequence 

dxu' -^ dxU weakly-* in [L°^((0, T); LlogL(J))]^ . (6.37) 

Form ( 16.36P and (I6.37p . we see that we can pass to the limit in the non-linear term in 
the sense 



[L\{Q,T)-Eexp{I))\'' -strong x [L°^{{{],T)-L\ogL{I))f-weak-^. 
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Because this is true for any bounded open interval / and for any T > 0, we deduce that, 

a{u') o d^u" -^ a{u) o d^u in V'{{0, T) x M) 
Consequently, we can pass to the limit in ^P^ and get that. 



dtu + a{u) o d^u = in ©'((0, T) x R). 

This solution u is also satisfy the following estimates (see for instance Brezis [3, Prop. 

3.12]): 

{El) ||f?xM||[ioo((o^T);LlogL(]R))]^ — -^™ ™ ll^i'""^!! [L°°((0,T);LlogL(IR))]" — ^' 
[E2) ||M||[ioo((o^T)xM)]*^ — ■^™i^MP'^ll[L°°({0,T)xR)]^^ — lpo|| [ioo(]R)]M , 

At this stage we remark that, thanks to these two estimates we obtain that {a{u) o 
dxu) e [L°°((0,T); L log L(M))] , which gives, since dtU = —a{u)odxU, that dtU G 
[L^ {{0, T);L log L{R))f\ and then u E [C([0,T); LlogL(M))]*^ 

Step 2 (The initial conditions): 

It remains to prove that the initial conditions ( IIP! ) coincides with m(-,0). Indeed, by 
Corollary l5.5l we see that, for all open bounded interval I of M, u^ is e-uniformly bounded 
in 



<M 



[w''\{o,Ty,w-''\i))Y' ^ cH[o,Ty,w~''\i)) 



M 



where W ^'^(/) is the dual of W^'°°{I). It follows that, there exists a constant C 
independent on e, such that, for all t,s E [0,T): 

\\u%t) - m"(s)||[^_i,i(^)]M < C\t - s|i 

In particular if we set s = 0, we have: 

h'W-«oll[VF-i.i(7)]*^<Cti (6.38) 

Now we pass to the limit in (16.381) . Indeed, the functions m"^ and Uq are e-uniformly 
bounded in [W^''^{{0,T);W~^'^{I))] and [H^^^'^(/)] respectively. Moreover we know 
that u'^ — Uq converges weakly-^ir in [L°°((0,T) x /)] to m — uq- 

Therefore, we can extract a subsequence still denoted by u'^ — Uq, that weakly-T*r converges 
in [Vr^'^((0, T); W~^'^{I))] to u — Uq. In particular this subsequence converges, for all 
t G (0,T), weakly-* in [L°°{{0,t);W~^'^{I))] , and consequently it verifies (see for 
instance Brezis 13, Prop. 3.12]), 
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\u - Mo||[ioo((o,i);vy-i,i(/))]A-f < liminf ||m^ - Uo\\[L^{{o,ty,w-^''^{i))]'^^ - ^'^'^ ■ 



From (16.380 we deduce that 

\\u{t) -Mo||[w'~i.i(7)]A^ < Ct 

which proves that m(-, 0) = uq in [P' 
Step 3 (Continuity of solution): 



M 



Now, we are going to prove the continuity estimate (II .Sp . For all /i > and (t,a;) G 
(0,T) X M, we have: 



[^(t,^ + /l) — M(t, x)| < 



x+h 



d^u{t,y)dy 



^ 2|| 1 II EXP(s,x+/i) il f^x'W|| Llog L(x,a;+/i) , 



< 2 



1 



<C 



Hi + 1) 
1 



|<9xm| 



L°°((0,T);LlogL(K)), 



Mi + 1)' 



where we have used in the second line the generalized Holder inequality (see Lemma 
16.41) and in last line we have used that dxU G L°°((0,T); LlogL(M)). Which proves 
finally the continuity in space. Now, we prove the continuity in time, for all (5 > and 

{t,x) G (0,T) X M, we have: 



5\u{t + 6,x) — u{t,x)\ 



x+5 



\u(t + 6,x) — u(t, x)\dy, 



< 



' l-x+5 

j u{t + 6,x)- 

K2 


u{t + S,y)\dy 


rx+S 

+ / \u{t + 6,y) 

Jx 


-u{t,y)\dy, 



+ 



x+5 



\u{t,y) -u{t,x)\dy. 



Similarly, as in the last estimate, we can show that: 
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x+5 fX+S 



px+o PX+O 

Ki + Ks <6 \d^u{t + 5,y)\dy,+5 1 \d,u{t,y)\dy. 

J X J X 

< 45||l||£;XP(x,x+<5)||t^xM|U°=({0,T);LlogL(IR)), 

-^Mi + i)- 

Now, we use that u is a solution of (JP]), and we obtain that: 



r-x+6 ft+b 

K2 < / \dMs,y)\dy, 

Jx Jt 



/t+5 px+5 
/ \a{u{s,y))od^u{s,y)\dsdy, 

< SMo\\u\\L°°{{0,T)xR)\\M\EXP{x,x+5)\\dxU\\L^((o,Ty,LlogL{M.), 

-^MJ + D- 

where we have used in last line that u G L°°((0, T^^^^ collecting the estimates of i^i, 
K2 and i^3, we prove that: 

\u{t^b,x)-u{t,x)\ <-{K, + K2 + Ks)<C , 

(5 ln(^ + 1) 

which proves finally the following: 

|M(t + 6,x + h) - u{t, x)\<C 



ln(i + l) ln(i + i: 



D 



7 Some remarks on the uniqueness 

In this Section we study the uniqueness of solution of the system (lPll- (lIDIl with 

j=l,...,M 

We show some uniqueness results for some particular matrices with M > 2. 
For the proof of Theorem ll.Sl in Subsection l7.2l we need to recall in the following Subsec- 
tion the definition of viscosity solution and some well-known results in this framework. 
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7.1 Some useful results for viscosity solutions 

The notion of viscosity solutions is quite recente. This concept has been introduced by 

Crandall and Lions [IHl [H] in 1980, to solve the first-order Hamilton-Jacobi equations. 

The theory then extended to the second order equations by the work of Jensen [27] 

and Ishii [23]. For good introduction of this theory, we refer to Barles [5] and Bardi, 

Capuzzo-Dolcetta [3]. 

Now, we recall the definition of viscosity solution for the following problem for all < 

e<l: 

dtV + H{t,x,v,d^v)-ed^^v = with x,v eR, t e {0,T). (7.39) 

where H : (0, T) x M^ i — > M is the Hamiltonian and is supposed continuous. We will set 

USC{{0,T) X M) = {/ such that / is upper semicontinuous on (0,T) x M}, 
LSC{{0,T) X M) = {/ such that / is lower semicontinuous on (0,T) x M}. 

Definition 7.1 (Viscosity subsolution, supersolution and solution) 

A function v G USC{{0,T) x M) is a viscosity subsolution of 1^7. 3 9\) if it satisfies, for 
every (to, a^o) ^ (0, T) x M and for every test function G C^((0, T) x M), that is tangent 
from above to v at {tQ^Xo), the following holds: 

dt(t> + H{to, xo, V, d^(t>) - ed^^cj) < 0. 

A function v G LSC{{Q, T) xW) is a viscosity supersolution of ^7.39^) if it satisfies, for 
every (to, xq) G (0, T) xR and for every test function G C^((0, T) x M), that is tangent 
from below to v at {to,XQ), the following holds: 

dt(f) + H{to, Xo, V, d^(f)) - ed^^cj) > 0. 

A function v is a viscosity solution of ^7.39\) if, and only if, it is a sub and a supersolution 

offjM). 

Let us now recall some well-known results. 

Remark 7.2 (Classical solution-viscosity solution) 

If V is a C"^ solution of ^7.39\) . then v is a viscosity solution of Il7.39i) . 

Lemma 7.3 (Stability result, see Barles |3, Th 2.3]) 

We suppose that, for e > 0, v^ is a viscosity solution of l(7.39i) . If v'^ ^ v uniformly on 

every compact set then v is a viscosity solution of ^7.39\) with e = 0. 

Lemma 7.4 (Gronwall for viscosity solution) 

Let V, a locally bounded USC{0, T) function, which is a viscosity subsolution of the 

equation —v = av where a > 0. Assume that v{0) < Vq then v <Vq e""^ in (0, T). 
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The proof of this Lemma is a direct application of the comparison principle, (see Barles 
[HI Th 2.4]). 

Remark 7.5 

From Lemmata \7.S\, 17.31 and from I6.35\) . we can notice that the solution m* of our system 
( f^j given in Theorem \l.l\ is also a viscosity solution of (^ (where the u^ for j ^ i are 
considered fixed to apply Definition \7.1\) . 

7.2 Uniqueness results 

In this Subsection we prove Theorem 11.51 Before going on, we recall in the following 
Remark a well-known uniqueness results and we recall in Theorem 17.71 the uniqueness 
results of W^'^ solution of (IF]). 



Remark 7.6 (Uniqueness for quasi-monotone Hamiltonians) 

// the elements of the matrix A satisfy: 

An + ^ A,j > for all z = 1, ■ ■ ■ , M. 

and if dxU^ > for i = 1, . . . , M, then we can easily check that the Hamiltonian 

Hi{u,d:,u')={ Y^ Aijund:,u\ 

\i=l,...M / 

is quasi-monotone in the sense of Ishii, Koike f2R (A. 3) J. Then the result of Ishii, Koike 
f2R Th.4-7j shows that for any initial condition uq G [L°°(]R)]*-'^ satisfying {H1)-{H2), 
the system (Hj satisfies the comparison principle which implies the uniqueness of the 
solution. 

We have the following result which seems quite standard: 

Theorem 7.7 (Uniqueness of the W^'°° solution) 

Let Mo G [Vr^'°°(M)]^^ and T > 0. Then system ^-^^) admits a unique solution in 

[W^'^{[0,T) xR)f. 

The proof of this Theorem is given in Appendix, because we have not found any proof 
of such a result in the literature. 

Proof of Theorem 11.5b 

Using Theorem 17.71 with a\u) = >^ Aiju\ it is enough to show that the system (JP))- 

j=l,...,M 

(llDll admits a solution in [iy^'°°([0,T) x R)] . To do that, it is enough to prove that 



the solution m"^ of the approximated system obtained in Corollary 15. 51 satisfies that d^u^ 
is bounded in [L°°((0,T) x M)] uniformly in < e < 1. Indeed, we then get the same 
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property for dxU, where u is the limit oi u^ as e ^ 0. Moreover, from the equation ((F 
satisfied by u and the fact that 

u e [L°°((0, T) X R)]^ and d,u E [L°°((0, T) x R)]^ , 
we deduce that dtu G [L°°((0,T) x R)]^^ which shows that u G [iy^'°°([0,T) x R)]^ . 

To simplify, we denote 

w'^ = dxU^-, 



and we interest in the 

maxw7^'*(t, x) = niiit). 



x& 



This maximum is reached at least at some point a;j(t), because tt;^'* G C°°((0,T) x R) fi 
Tyi'P((0,T) X R) for all 1 < p < +00 (see Lemma^H Km ). 

In the following we prove in the two cases (i) and (ii) defined in Theorem 11.51 that 
mj, for alH = 1, . . . , M, is bounded uniformly in e. First, deriving with respect to x 



the equation ( [^]) satisfied by u^ G [C°°((0,T) x R)] , we can see that w^ satisfies the 



following equation 

dtw''' - edxxw''' + Yl Aju''^dxw'''+ ^ AiX'^w''' = 0. (7.40) 

j=l,...,M j=l,-,M 

Now, we prove that nii is a viscosity subsolution of the following equation, 

jmi{t)+ J2 Aijw'^^{t,Xi{t))w''\t,Xi{t))<0. (7.41) 

i=l,...,M 

Indeed, let cj) G C^(0,T) a test function, such that cj) > rui and (^(to) = "^■i(^o) for some 
to G (0,T). From the definition of mj, we can easily check that > w^'''{t,x) and 
0(to) = w^'\to,Xi{to)). But, the fact that w''' G C°°((0,T) x R), by Remark O we 
know that w^'^ is a viscosity subsolution of ( 17.40p . We apply Definition 17. 11 and the fact 
that dx(p = dxx<P = 0, we get 

-0(to)+ Yl AX''(^o,a;.(to))«^''Xto,a;i(to))<0. 

j=l,...,M 

Which proves that rrii is a viscosity subsolution of (17.411) . 

Two cases may accur: 

i) Here, we consider the case where M > 2 and Aij > for all j > i. We see the 
equation satisfied by mi, we deduce that satisfies (a viscosity subsolution) 
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j=l,...,M 

where we have used the fact that, for j = 1, . . . , M, Aij > and w'^'^ > 0. This proves 
by Lemma EH (with a = 0) that, 

mi(t) < mi(0) = w^'^(t,xi(t)) < ||9^Mo||loo(r). 

We reason by recurrence: we assume that rrij < C for all j < i, where C is a positive 
constant independent of e, and we prove that mj+i is bounded uniformly in e. Indeed, 
we know that 

-m,+i(t) <- Yl A+i,jw''\t,Xj{t))w'^'^\t,Xi+iit)), 

j=l,...,M 

<- J2 A+i,jW'''it,x,it))w'''+\t,x,^iit)) 

- Y, A+l,jW'''{t,Xj{t))w'''+\t,Xi+i{t)), 
M>j>i+1 

We use that ^j+ij > 0, for M > j > i + 1, we obtain that 



<C V \Ai+ij\ mi+i(t). 






This implies by Lemma [7\4l with a = C y^ \Ai^ij\ ) , that 

Which proves that for alH = 1, . . . , M, rrii is bounded uniformly in e. 

ii) Here, we consider the case where M > 2 and Aij < for all i j^ j. Taking the sum 
over the index i, from (I7.4ip we get that the quantity m{t) = > "^j(t) satisfies (a 



=1,...,M 



viscosity subsolution see Bardi et al. |^) 
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i,j=l,...,M 

<- ^ AiX'^{t,x,{t))w''\t,Xi{t)), 

i,j=l,...,M 

<0. 

where we have used that the matrix A satisfies {H2') and w^'* > 0, for i = 1, . . . ,M. 
Using Lemma EH with a = 0, we get 



i{t) <m(0)= Yl ^- 



m[t] <^ m(U) = > OxUq\ 



i=l,...,M 

< sup V d^UQiy). 

y^^i=i,...,M 

which proves l\1.7\ i. □ 

8 Application on the dynamics of dislocations densi- 
ties 

In this Section, we present a model describing the dynamics of dislocations densities. 
We refer to [22] for a physical presentation of dislocations which are (moving) defects 
in crystals. Even if the problem is naturally a three-dimensional problem, we will first 
assume that the geometry of the problem is invariant by translations in the Xa-direction. 
This reduces the problem to the study of dislocations densities defined on the plane 
(a;i,X2) and propagation in a given direction b belonging to the plane (xi,X2) (which is 
called the "Burger's vector"). 

In this setting we consider a finite number of slip directions b E M.'^ and to each b we 
will associate a dislocation density. For a detailed physical presentation of a model with 
multi-slip directions, we refer to Yefimov, Van der Giessen [Hj and Yefimov [401 ch. 5.] 
and to Groma, Balogh [2T] for the case of a model with a single slip direction . See 
also Cannone et al. [8] for a mathematical analysis of the Groma, Balogh model. In 
Subsection ??, we present the 2D-model with multi-slip directions. 

In the particular geometry where the dislocations densities only depend on the variable 
a; = xi + X2, this two-dimensional model reduces to one-dimensional model which pre- 
sented in In Subsection I8.2[ See El Hajj [15] and El Hajj, Forcadel [16] for a study in 
the special case of a single slip direction. Finally in Subsection 18.31 we explain how to 
recover equation (JP]) as a model for dislocation dynamics with a\u) = N^ AijU^ for 

j=l,...,M 

some particular non-negative and symmetric matrix A. 
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8.1 The 2D-model 

We now present in details the two-dimensional model. We denote by X the vector 
X = (xi,X2). We consider a crystal filling the whole space M? and its displacement 
V = (fi,f2) : K^ — ^ IR^, where we have not yet introduced the time dependence for the 
moment. 

We define the total strain by 

e{v) = -iVv + ^Vv), 

where Vv is the gradient with (Vv)ij = — ^, i,j G {1, 2}. 

OXj 

Now, we assume that the dislocations densities under consideration are associated to 
edge dislocations. This means that we consider M slip directions where each direction 
is caraterize by a Burgers vectors b'' = {b\, 62) G IR^, for /c = 1, . . . , M. This leads to M 
type of dislocations which propagate in the plan (xi,a;2) following the direction of b'^, 
iov k = l,...,M. 

The total strain can be splitted in two parts: 

£{v) = e' + e". 
Here, £^ is the elastic strain and e^ the plastic strain defined by 



e^ 



J2 ^°''^w^ (8-42) 

k=l,...,M 



where, for each k = 1, . . . , M, the scalar function u^ is the plastic displacement associated 
to the k-th slip system whose matrix e^'^ is defined by 



0,fe 

2 



^(V'^n'' + n^^b'y (8.43) 

where n*^ is unit vector orthogonal to r and Ib'^ (^ ft''] = b\i 



^3 



To simplify the presentation, we assume the simplest possible periodicity property of 
the unknowns. 

Assumption (H): 

i) The function v is Z"^ -periodic with / v d'X. = 0. 

i(0,l)2 

a) For each k = 1, . . . , M, there exists L^ G M^ such that u^ — L^ -IL is a 1? -periodic. 
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Hi) The integer M is even with M = 2N and L^+^ = L^, and that 



'k+N 



L\ P+ 



N 



-.0,k+N 



_0,A: 



n 



e ' = —e 

iv) We denote by f'' = {ti,T2) a vector parallel to h^ such that r'^+^ 
that L^ is chosen such f^ ■ L'^ > 0. 



f^. We require 



The plastic displacement u^ is related to the dislocation density associated to the Burgers 
vector h^ . We have 



/c-th dislocation density = r^ ■ Vu^ > 0. 

The stress is then given by 

o- = A : e^ 

i.e. the coefficients of the matrix a are: 



(8.44) 
(8.45) 



a 



u 



^ Kijkieli for ?,j = 1,2, 



k,l=l,2 



where A = {^ijki) i j ^ i=i 2^ ^^e the constant elastic coefficients of the material, satisfying 
for m > 0: 

^ Ai,j,k,i£ij£ki > m ^ efj (8.46) 

ijkl=l,2 i,j=l,2 

for all symmetric matrices e = {£ij),j, i.e. such that eij = e 



^])ij^ 



-y "-J*- 



Finally, for fc = 1, . . . , M, the functions u^ and v are then assumed to depend on 
(t,X) G (0,T) X M^ and to be solutions of the coupled system (see Yefimov [40l ch. 5.] 
and Yefimov, Van der Giessen [4T]): 



div a 

a 

e{v) 





A : {e{v) - e'P) 
\ {Vv + ^Vv) 

fc=l,...,M 



on (0,r) xM^ 
on (0,T) X M^ 
on (0,T) X M^ 
on (0,T) xR^, 



^.47) 



(^ dtu^ = {a : e^^^)T^.Vu^ on (0, T) x M^ for fc = 1, . . . , M, 
i.e. in coordinates 
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^ dxi 

7 = 1,2 J 



0", 



^ Aijkiiskiiv) 



e,. 



k,l=l,2 

, , 1 / dvi dvj ~ 

0,fc„,A: 



on(0,T)xM^ fori = 1,2, 
on (0,T) X M^ 



£: 



E 

fc=l....,M 



£,, n 




an'^ 



on (0,T) X R^ 
on (0,T) X M^ 

t^Vm'^ on(0,T)xl 



for z,j = 1,2 



(8.48) 



for A; = 1, . . . , M, 



where the unknowns of the system are u^ and the displacement v = (fi,f2) and with 
e^'^ defined in (18.431) . Here the first equation of (18.471) is the equation of elasticity, while 
the last equation of (18.471) is the transport equation satisfied by the plastic displacement 
whose velocity is given by the Peach- Koehler force cr : e^''^. Remark that this implies 
in particular that each dislocation density satisfies a conservation law (see the equation 
obtained by derivation, using (18. 44^ ). Remark also that our equations are compatible 
with our periodicity assumptions (H), {i)-{ii). 



8.2 Derivation of the ID-model 

In this Subsection we are interested in a particular geometry where the dislocations 
densities depend only on the variable a; = xi + a;2. This will lead to ID-model. More 
precisely, we make the following: 

Assumption {H'): 

i) The functions v{t, X) and u'^{t, X) — L*^ ■ X depend on the variable x = X\^ X2- 

%i) t\ + r|' = 1, /or A; = 1, . . . , M. 

Hi) L\ = L\jork = \,...,M. 



For this particular one-dimensional geometry, we denote by an abuse of notation the 
function v = v{t,x) which is 1-periodic in x. If we set /* 



''^ — ^^ % we have 



L^-X = l'' ■x + 



J k J k' 

^ '^ ^ {Xi -X2). 
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By assumption {H'), {m), we see (again by an abuse of notation) that u = {u^{t, x))k=i,...,M 
is such that for A; = 1, ... , M, u^(t, x) — l^ ■ x is 1-periodic in x. 

Now, we can integrate the equations of elasticity, i.e. the first equation of (18 .47^ . Using 
the Z^-periodicity of the unknowns (see assumption (if), (z)-(m)), and the fact that 
^o,k+N _ _^o,k ^ggg assumption (if), (m)), we can easily conclude that the strain 



e^ as a linear function of {u^ — M-'^^)j=i,...,Ar and of I / {u^ — u^^^) dx J 

This leads to the following Lemma 

Lemma 8.1 (Stress for the ID-model) 

Under assumptions (H), {i)-{ii)-{iii) and {H'), {i)-{iii) and ([g.^gp , we have 



j=l,...,N 

(8.49) 



y^ AjU^ + ^ Qij I u^ dx, fo 

=1,...,M 7=1,. ..,M ° 



- a : £°'^ = > J Aiju^ + > J Qij / u^ dx, for i = I, . . . ,N. (8.50) 

j= 

where for i,j = l,...,N 

(8.51) 

Qi,j = Qj,i C'lT'd Qi+NJ = ~Qi,j ~ Qi,j+N- 

Moreover the matrix A is non-negative. 

The proof of Lemma [8?T] will be given at the end of this Subsection. 

Finally using Lemma ISTTl we can eliminate the stress and reduce the problem to a one- 
dimensional system of M transport equations only depending on the function u\ for 
i = 1, . . . , M. Naturally, from (I8.50p and (if') 7 (^0 this ID-model has the following form 



dtu^+ y^ AijU^ + Y^ Qij / u^ dx I drcU^ 
\j=l,...,M j=i,-,M >^o / 



with from (I8.44p 



0, on(0,T)xM, for i = 1,...,M, 

(8.52) 



d.^u'>0 for z = 1,...,M. (8.53) 

Now, we give the proof of Lemma ISTTl 
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Proof of Lemma 18.1b 

For the 2D-model, let us consider the elastic energy on the periodic cell (using the fact 

that e^ is Z^-periodic) 

E^^ = - [ A:e':e^ dX. 
2 i(o,i)2 

By definition of a and e^, we have for i = 1, . . . , M 

a : e°'^ = -V„,E"'. (8.54) 

On the other hand usind {H'), {i)-{iii), (with x = Xi + X2) we can check that we can 
rewrite the elastic energy as 

1 r^ 

E"^ = - K:e^:e^ dx. 



2./0 



Replacing e^ by its expression (18.491) . we get: 



•'^ i,i=h...,N 



E^^= - I y Aij{u^ - u^+^){u' - u'+^) dx 






;,j=i,...,iV ^- " . .- 

for some symmetric matrices Aij = Aj^i, Qij = Qj^i. In particular, joint to (I8.54p this 
gives exactly ( 18301) with (EM]). 

Let us now consider the functions w* = u* — m*+^ such that 

-1 

w^ rfx = for i=l,. . . ,N, (8.55) 

'0 

From (18.46P that we deduce that 



1 r 

0<E'^ = - V Aijw'w^ dx. 

2 Jo .. , 1 



ij=l,...,Af 

More precisely, for alH = 1, . . . ,N and for all w^ G M, we set 

^^^iw' on [0,i], 

\ -w' on [|,1], 

which satisfies (I8.55p . Finally, we obtain that 



u: E A 



< E-^' = - / > A.w^w^ rfx. 



2 / ^-^ *"' 

"^0 ij=l,...,Af 

Because this is true for every w\ we deduce that A a non-negative matrix. D 

37 



8.3 Heuristic derivation of the non-periodic model 

Starting from the model (18.52I) - (I8.53I) where iov i = 1, . . . , M„ u''{t, x) — /*-a; is 1-periodic 
in X, we now want to rescale the unknowns to make the periodicity disappear. More 
precisely, we have the following Lemma: 

Lemma 8.2 (Non-periodic model) 

Let u be a solution of Ii8.52) - ^.5S\} assuming Lemma lKll and u\t, x) — l^ -x is 1-periodic 
in X. Let 

u-g(t, x) = u^{6t, 6x), for a small 6 > and for j = 1, . . . , M, 
such that, for all j = 1, . . . , M 



ui{0,-)^u^{0,-), as 6^0, and u\0,±oo) = u^^^{0,±oo) (8.56) 

Then u = {u-^)j=i,...,m formally is a solution of 

dtu'+i Yl AjuAd^u' = 0, on(0,T)xM, (8.57) 

\j=l,...,M ) 

with the matrix A is non-negative and d^u^ > for i = 1, . . . , M. 

We remark that the limit problem (I8.57P is of type ([P]) with (i^l') and {H2'). 



Now, we give a formal proof of Lemma [8l 

Formal proof of Lemma 18. 2t 

Here, we know that u\ — 5r ■ x is --periodic in x, and satisfies for i = 1, . . . , M 



dt< + ( 5^ A A + ^ 5Z Q'l I' "5 ^^ ) ^-"5 = 0, on (0, T) X M, (8.58) 

\i=l,...,M J=1,...,M -^0 / 

To simplify, assume that the initial data ^^(O, ■) converge to a function m(0, ■) such that 
dxUs[Q,-) has a support in {—R,R), uniformly in 6, where R a positve constant. We 
expect heuristically that the velocity in (I8.58P remains uniformly bounded as 5 — *> 0. 

Therefore, using the finite propagation speed, we see that, there exists a constant C 
independent in 6, such that dxUs{t, ■) has a support in {—R — Ct, R + Ct) uniformly in 
5. Moreover, from f l8.56p and the fact that 

y^ Qij / u{ dx = ^ Qij / {u^ - u^^^) dx, 

j=l,...,M "^0 j=l,...,N "^0 
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we deduce that 



Y^ Qij / ul dx, 

-1 A# JO 



''5 "'■^? 
j=l,...,M 

remains bounded uniformly in 6. Then formally the non-local term vanishes and we get 
for z = 1,...,M 



22 ^ij'^i + ^ z2 '5*-'' / "^i dx ^ 2J ^ijU^ , as 5 — > 0, 

i=l,..,M i=l,.-,M "^° J=1,-,M 

which proves that m is solution of (18.571) . with the matrix A is non-negative . □ 



9 Appendix: proof of Theorem 17.71 

Let Ml = {u])i and M2 = (^2)1, for i = 1, ■ ■ ■ , M, be two solutions of the system (lP|) in 
[Vr^'°°((0,T) xM)]*'^ such that Mi(0,-) =m^(0,-)- 

Then by definition u\ and U2 satisfy respectively the following system, for i = 1, ■ ■ ■ , M: 

dtu\ = -a'{ui)d^u\, 

Subtracting the two equations we get: 

dt {u\ - 4) = - {a'iui) - a\u2)) d^u\ - a\u2)d,^iu\ - 4)- 

Multiplying this system by {u\ — u^) {ipY where ipi^) = e"'^', and integrating in space, 
we deduce that: 

-17 p 



- / a\u2)^l)'^ ("i - 4) 5^(< - 4)- 



Taking the sum over i, we get: 



/i 



1 d 
2di 



:1,...,M / '^''^ J = 1,...,M 



^i=l,...,M 



/2 



'^''^i=l,...,M 



4)^ . 



39 



Integrating I2 by part, we obtain: 



hi 



■^^ i,i=l,...M 



I22 



•^^i=l,...,M 



U^fd.iij^) . 



Next, using the fact that M2 is bounded in H^^'°°((0, T) x M), for i = 1, . . . , M, we deduce 
that: 



I/21I < |MMi||M2||[^^((o,r)xKrH X] ||«-4)^| 

\J=1,...,M / 

(9.59) 



|2 
=1,...,M 



Since (9j:(-?/'(x))^ = —2sign{x){%l){x)Y and M2 is bounded in W^^'°°((0,T) x M), for i 
1, ■ ■ ■ , M, we obtain: 



L, )-"-') 



1^22! <\mA J2 



{u{-ui)iijf 



L2(I 
,...,M 

(9.60) 



<c E llM-«2)v^l 



|2 

Il2( 
vJ=l,...,M 



Now, using the fact that u\ is bounded in W^'°°{{0,T) x M), for i = 1,-, -jM, and the 
inequality \ab\ < |(a^ + 6^), we get: 



|/l| < -Mi(M+1)||mi||[h..o((o,t)xR)]M / Yl 1^1 -^2lV, 

'^'*i=l,...,M 



2 



< -Mi(M+ l)||Mi||[t^^((o,T)xM)]M E ||« -M*2)V^||l2(k) ' (9.61) 



vi=l,...,M 



<c E IIK-«2)^I 



|2 

Il2( 
U=1,...,M 
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Finally, flMTll . (19391) and (iOnH . imply: 

^^l ^ ||K-«*2)^||'.(^J <2(|Ji| + |J2i| + |/22|)<C7( 5^ ||K-^)^||'.(, 
\i=l,...,M / \i=l,...,M 

Now, we apply the Gronwall Lemma and we use that u\{0, ■) = ^(0' ')? to deduce that: 

J2 ||K-4)V^|ll.c((o,T);L^(M))< E IIH(0'-)-4(0,-))^||l.(R)e^^ = 0, 

1=1,. ..,M i=l,...,M 

i.e., Ml = U2 a.e in (0,T) x M. D 
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